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Abstract. We study a smooth cosmological solution of the 5D anisotropic standing 
wave braneworld model generated by gravity coupled to a phantom-like scalar field. In 
this model the brane emits anisotropic waves into the bulk with different amplitudes 
along different spatial dimensions. We found a natural mechanism which isotropizes 
the braneworld, rendering a 3-brane with de Sitter symmetry embedded in a 5D de 
Sitter space-time for a wide class of initial conditions. The resulting thick geometrical 
braneworld (a de Sitter 3-brane) possesses a series of remarkable features. By explicitly 
solving the bulk field equations we are able to give a physical interpretation of the 
anisotropic dissipation: as the anisotropic energy on the 3-brane rapidly leaks into the 
bulk, through the nontrivial Weyl tensor components, the bulk becomes less isotropic. 
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1. Introduction 

Braneworld models involving large extra dimensions have been very useful in addressing 
several open questions in high energy physics [U[2] , and astrophysics and cosmology (for 
reviews see [HHH]). For more than a decade it is clear that braneworlds possess a set 
of appealing features that encourage one to develop this research line further. Namely, 
braneworld models recast the hierarchy problem into a higher dimensional viewpoint, 
reveal a geometrical mechanism of dimensional reduction supported by a curved extra 
dimension, can trap various matter fields on the brane (including gauge bosons in 
some cases), take into account the self-gravity on the brane through its tension, give a 
realization of the AdS/CFT correspondence to lowest order, allow one to address the 
smallness of the cosmo logical constant, provide several cosmological backgrounds with 
consistent dynamics that lead to General Relativity results under suitable restrictions on 
their parameters, etc. Most of the branworlds were realized as time independent field 
configurations. However, mostly within the framework of cosmological studies, there 
have been proposed different braneworld models that employ time- dependent metrics 
and matter fields (if any), as well as branes with tensions varying in time [6til4]. 

It is widely believed that our Universe started its evolution from a highly anisotropic 
state which lead to a highly symmetric state that we observe nowadays (according to 
the recent WMAP data CMB is isotropic with an accuracy of 10~ 5 [IS]). An interesting 
cosmological question concerns the dynamical mechanism that led the Universe to shed 
almost all of its anisotropy along time till the present epoch which is highly isotropic. 
A very popular mechanism among cosmologists is inflation: the Universe diluted away 
almost all of its anisotropy during a period of accelerating expansion in an early epoch. 
Of course, there are several other mechanisms that tried to explain this important 
phenomenon (see [12] and references therein). 

When studying braneworlds generated by anisotropic backgrounds, one must have 
a complete solution to the bulk and brane field equations in order to consistently analyze 
the cosmological dynamics on the brane. Such a solution involves the knowledge of the 
bulk Weyl tensor which is felt on the brane by means of its projection [16]. It turns out 
that this is not an easy task to solve [9] and up to now, there are a very few cosmological 
anisotropic braneworlds that present a complete solution. In this context, there are 
many studies reported in the literature that made several assumptions about the Weyl 
term in the absence of knowledge of the full solution to the bulk metric p21[TE]. The 
main difficulty here consists in finding anisotropic generalizations of the AdS§ space-time 
that are necessarily non-conformally flat and incorporate anisotropy on the cosmological 
brane. In this sense, it is relevant to propose new braneworld generalizations that 
attempt to describe more realistic cosmological braneworld models, or explore other 
aspects of higher- dimensional gravity which were not addressed till now by the known 
brane models. 

Recently a braneworld model generated by 5D anisotropic standing gravitational 
waves minimally coupled to a phantom-like scalar field in the bulk was proposed [T9||20] . 
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In this framework, an alternative mechanism for localizing 4D gravity as well as for 
trapping matter fields [21], including gauge bosons [22], which usually are not localized 
on thin braneworlds, was found. 

In this paper we consider a thick version of the braneworld model [T91I20] and 
slightly modify the metric ansatz by adding a scale factor that multiplies the spatial 
sector of the anisotropic 3-brane in the 5D metric. When exactly solving the bulk 
field equations the cosmological constant turns out to be positive, realizing a de Sitter 
braneworld instead of an Anti de Sitter one. By leaving the found solution evolve in 
time, the metric isotropizes and the corresponding scalar field disappears, leading to a 
recently proposed de Sitter braneworld model purely generated by curvature, i.e., by an 
interplay between the 4D and 5D cosmological constants [23]. It turns out that within 
this regular braneworld one is able to model the thick brane in a completely geometrical 
way, avoiding at all the use of scalar matter. In this model the effective 4D Planck mass 
is finite (and also depends on the Hubble parameter of the metric), the 4D gravity can be 
localized and the correction to Newton's law have been computed. As an extra bonus, a 
mass gap is displayed in the gravity spectrum of Kaluza-Klein (KK) excitations (a fact 
that fixes the energy scale at which these massive fluctuations can be excited and enables 
us to avoid difficulties when analyzing the traces of ultra-light KK excitations) without 
developing naked singularities as in the case of scalar thick brane configurations |24j. 
Another interesting feature of this thick braneworld model is that its 4D cosmological 
constant can be made as small as one desires without the need of fine-tuning it with the 
bulk cosmological constant, as it happens in the Randall-Sundrum type models [2]. 



2. The model and the complete solution 

We start with a 5D action which describes gravity coupled to a non-self-interacting 
scalar phantom-like field [25], which depends on time and propagates in the bulk: 



S b = J d b Xy/g 



1 ,-..2 



(R - 2A) + - (V0) 



(1) 



_16ttG v 2 

where G and A are 5D Newton and cosmological constants, respectively. To avoid the 
well-known problems of stability which occur with ghost fields we assume that the bulk 
phantom-like scalar field does not couple to ordinary matter in the model. 
The Einstein equations for the action ([T]) are written in the form: 

12 2 
R^v = - -g^T + -g^A = -d^d v a + -g^A , (2) 



where Greek indices refers to 5D and we have redefined the scalar field as a = vo7rG0 
in order to absorb the 5D gravitational constant. 
Here we shall use the anisotropic metric ansatz: 

ds 2 = e 2A{w) [-dt 2 + a 2 (t) (e u dx 2 + e u dy 2 + e~ 2u dz 2 ) + dw 2 ] , (3) 

where a(t) plays the role of a cosmological scale factor, u(t,w) describes standing 
gravitational waves of [19,20] an d the warp factor A(w) is an arbitrary function of 
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the extra coordinate w, allowing for smooth solutions, and models a thick brane 
configuration in the spirit of [26J. This metric generalizes straightforwardly the thin 
brane metric ansatz of [T91I20] . to the case when A(w) ^ \w\ and a scale factor a(t) 
multiplies the 3D anisotropic spatial sector of the metric. 

Thus, in this braneworld, the 3-brane possesses anisotropic oscillations which send 
a wave into the bulk (as in [27]), i.e. the brane is warped along the spatial coordinates 
x,y,z through the factor e u ^' w \ depending on time t and the extra coordinate r. 
Moreover, these warped spatial coordinates are multiplied by a scale factor a(t) that 
allows them to evolve in time (expanding or contracting) as an extra feature. This 
model also joins previously constructed anisotropic braneworlds [8HT3] which approached 
several cosmological issues like anisotropy dissipation during inflation [10], braneworld 
isotropization with the aid of magnetic fields [12] and localization of test particles 
[TBI l2~T| [22] . Moreover, as a general feature it has been established that anisotropic 
metrics on the brane cannot be supported by static bulks [TT|[T2]. 

On the other hand, the phantom-like scalar field a(t, w) obeys the Klein-Gordon 
equation on the background (j3J): 

' (V^gg^dvO-) = a + — & - a" - 3A'a' = , (4) 



V~9 a 
where overdots and primes stand for derivatives with respect to time and extra 
coordinates, respectively. 

Under the ansatz ([3]) we can write the Einstein equations in the form: 

Ru = -h 2 -3-+A" + 3A' 2 = --e 2A A-& 2 , 
2 a 3 



2 M 



- w + 3-«-/-3iV]+- + 2n - A "-3A' 2 
2 \ a J a \a 



2 



R yy = \e 2A+u a 2 k 



-u - 3-u + u" + 3A'u' + - + 2 [ - ] - A" - 3 A' 2 



2 



(5) 



R zz =a 2 e~ 2u 

= l a 2 e 2A-2u A 

3 

R ww = - \u ,2 -AA"= 2 -e 2A A-o- 12 , 
3 

Rtw = - ~uu' = -era 1 . 

From the tw-component of the Einstein equations (jSJ), it follows that the fields a 
and u are related by 

a(t ' r) = vf u(t ' r) ' (6) 

up to a meaningless additive constant (the same relationship was obtained in [T9|[20]). 
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We further take into account the Klein-Gordon equation (jl]) in the Einstein 
equations (jSJ), and combine the tt- and xx- components to obtain a single relevant 
equation for the scale factor: 

a \a J 
which has a solution of the form 

a{t) = a e m , (8) 

where a$ is a constant that can be absorbed in a redefinition of the spatial coordinates. 

Once we know the solution for the scale factor, we further substitute its expression 
back into the Einstein equations, and combine the xx- and ww-components in order to 
get an equation for the warp factor A(w): 

A" - A' 2 + H 2 = . (9) 

This equation has the solution: 

H 



A(w) = In 



-sech [H(w — wq)} 



(10) 



where the integration constant b is related to the 5D cosmological constant as follows: 
A 5 = 66 2 , (11) 

which reveals the de Sitter nature of the 5D bulk space-time. 

On the other hand, after taking into account the solution for the scale factor (jHJ), 
and assuming the following ansatz for the scalar field: 

a{t, w) ~ u(t, w) = e(t)xHe-i A , (12) 

the Klein-Gordon equation reduces to a couple of ordinary differential equations: 

X ~ + 9 -A' 2 -Q^ X = 0, (13) 

e + 3He + tt 2 e = , (14) 

where Q is an arbitrary constant. 

The equation ([TBI can be regarded as a Schrodinger equation with the analog 
quantum mechanics potential: 

Vqm = \a" + \A' 2 . (15) 



Using the solution for the warp factor ( JTOl) we see that it represents a modified Poschl- 
Teller potential of the form: 

V QM (w) = ^H 2 - ^tf 2 sech 2 [H(w - w )] . (16) 

Thus, the differential equation for x (O turns out to be a known eigenvalue problem 
with a mixed spectrum (see for instance). Namely, there are two bound states: a 
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ground state with Q = and another one with Q = \^2H, separated by a mass gap 
determined by the asymptotic value: 

V QM (oo) = hi 2 , (17) 

as well as a continuum of KK states starting at Q = 3H/2. 

In fact, the equation for x is precisely the eigenvalue equation for the transverse 
traceless KK metric fluctuations and massless scalar perturbations of the recently 
proposed de Sitter braneworld model [221121], where there is a mass gap in the graviton 
and scalar mass spectra of KK perturbations. 

Thus, the equation for the relevant function x possesses a general solution: 

X (w) = [tanh (H(w - w ))} + C 2 Q% /2 [tanh (H(w - w ))) , (18) 

where P^ 2 and Q^ 2 are associated Legendre functions of first and second kind, 
respectively, grade v = 3/2 and order 



9 n 2 

" = Vi-iP- (19) 

Thus, the first discrete state is the ground state with Q = 0, order /i = 3/2 and 

E = ~H 2 . (20) 

The explicit expression of this zero mode is given by: 

X o(w) = c sech 3/2 [H{w - w )} , (21) 

where Cq is a normalization constant. On the other hand, the second bound state 
corresponds to an excited mode with f2 = \/2H, order /i = 1/2 and energy 

Ex = -\h 2 , (22) 
and has the following form: 

Xi{w) = ct sinh [H(w - w )} sech 3/2 [H{w - w )} , (23) 

where C\ also is a normalization constant. 

Turning back to the equation for e(t) ( !T4l) we see that it is the equation for a damped 
oscillator which has three different solutions depending on the kind of damping, i.e. on 
the relation between the constants H and Q: a), under-damping (Q 2 > 9H 2 /4), b). 
critical damping (Q 2 = 9H 2 /4), and c). over-damping (Q 2 < 9H 2 /4). 

a). The solution for the under-damped case reads: 

e(t) = Ce~^ m sin {ut + 5) (24) 
where the constant C denotes the oscillations amplitude, the parameter 



uj = x /n 2 - 9 -H 2 (25) 

is the un-damped frequency of the oscillations and 5 is an arbitrary phase constant. 
The constants C and 5 are determined by initial conditions. We can easily see that 
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these oscillations will exponentially decay to zero with time, a fact that translates into 
a vanishing metric function u — Y 0, which leads in turn to an isotropic 5D metric with 
a 3-brane with de Sitter symmetry [23] . 

b) . The solution for the critical damping case reads: 

e(t) = e -'^ Ht (at + ft) , (26) 

where a and /3 are arbitrary constants determined by initial conditions. In this case we 
observe the same effect as in the under-damped one: the amplitude of the function e 
exponentially vanishes with time, making the metric function u disappear and yielding 
a de Sitter 3-brane embedded in an isotropic 5D de Sitter space-time. 

c) . Finally, the over-damped case possesses the following solution: 

e(t) = e~i Ht (ae M + Pe~ M ) , (27) 

where 

u = \J\h 2 ~W , (28) 

and a and are arbitrary constants determined by initial conditions. Once again we 
see that the metric function u exponentially decays in time, leading to a completely 
isotropic 5D metric with a de Sitter 3-brane embedded in it and realizing a very natural 
isotropization mechanism. 

Thus, we see that in the general case, the solution for the time evolution of the 
metric function u(t,w) expressed by ( fl2l exponentially yields an isotropic 5D metric 
which possesses an embedded 3-brane with de Sitter symmetry under a wide class of 
initial conditions. It is worth noticing that together with the function u, the scalar 
field also exponentially disappears as a consequence of their proportionality, rendering 
a completely geometric de Sitter thick brane with very appealing features (see [231 129] 
for details). 

In order to understand in detail what happens to the anisotropy of the initial metric 
and give a consistent physical interpretation to this process, it is necessary to compute 
the Weyl tensor and its projection to the brane using, for example, the formalism of [16J. 
However, it is clear that the inflationary process that isotropizes the metric is due to 
the 4D cosmological constant (since it is equal to A4 = 3H 2 ), which is also responsible 
for the vanishing of the phantom-like bulk scalar field. 

3. Conclusions 

In this paper we have presented a novel mechanism of isotropization of an initially 
anisotropic 5D thick braneworld generated by a phantom-like scalar field minimally 
coupled to gravity. Under a wide class of suitable initial conditions the anisotropic 
braneworld evolves and exponentially isotropizes by itself, under the action of the 4D 
cosmological constant, while the scalar field disappears. Thus, the anisotropic energy of 
the 3-brane rapidly leaks into the bulk through the nontrivial components of the Weyl 
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tensor, leading to a less isotropic bulk. As a result we have an isotropic thick braneworld 
supported by pure curvature, where 4D gravity as well as other matter fields have been 
shown to be localized, the effective 4D Planck mass is finite and depending on the 
Hubble parameter, the 4D cosmological constant can adopt an infinitely small value 
without any fine-tuning to the 5D one, and the mass spectrum of KK gravitons displays 
a mass gap, an important feature from the phenomenological point of view; moreover, 
the corresponding corrections to Newton's law have also been computed [23l[29] . 
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